ALGORITHMIC COOLING OF A QUANTUM SIMULATOR 
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Abstract. Controlled quantum mechanical devices provide a means of simulating more complex 
quantum systems exponentially faster than classical computers. Such "quantum simulators" rely 
heavily upon being able to prepare the ground state of Hamiltonians, whose properties can be used 
to calculate correlation functions or even the solution to certain classical computations. While 
adiabatic preparation remains the primary means of producing such ground states, here we provide 
a different avenue of preparation; cooling to the ground state via simulated dissipation. This is in 
direct analogy to contemporary efforts to realize generalized forms of simulated annealing in quantum 
systems. 
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1. Introduction. Quantum devices provide new opportunities in communica- 
tion and computation [1, 2, 3, 4, 5, 6, 7, 8, 9]. One promising application of a well 
controlled quantum device is simulating a quantum system, which can occur expo- 
nentially faster than can be achieved classically [2, 10, 11, 12]. Such simulations could 
provide insights in many current fields of research, such as BCS-BEC supcrfluids 
[13, 14], quantum chemistry [15, 16], and highly correlated condensed matter systems 
[17, 18, 19, 20]. However, a crucial component of such simulation is the specification 
of the initial state of the system to be simulated. While for some problems, including 
many in quantum chemistry [21], such initial states can be found via prior knowledge 
from classical computer studies, in general a means of preparing such states does not 
exist. 

Methods for the preparation of specific eigenstates of Hamiltonians, particularly 
the ground state, therefore remain a pressing challenge for the most interesting quan- 
tum simulation applications. One approach for preparation is adiabatic evolution 
from a system with an accessible ground state [22]. Here, we offer an alternative ap- 
proach: cooling to the ground state by expanding the simulation to include a (small) 
quantum bath. This approach differs crucially from prior work in several respects. In 
contrast with adiabatic approaches [23], it requires only information concerning the 
spectrum of the cooled Hamiltonian Hs and not any intermediate Hamiltonians of the 
form XHs + (1 — A)7?o- In its most general form it is able to prepare ground states of a 
wide class of gapped Hamiltonians. Specifically, it may cool any gapped Hamiltonian 
with a tight band of excited state energies that are separated from the rest of the 
spectrum. For appropriate cases, we find our approach provides a quadratic speedup 
for ground state preparation in analogy with Grovcr's algorithm [24]. Finally, the 
approach may also be used to cool a variant of Kitaev's clock Hamiltonian [25] at a 
polynomial overhead, so that it may efficiently solve any problem solveable through 
conventional circuit-based quantum computation [26] . 

2. Grover's Algorithm by Simulated Cooling. As an illustrative example, 
we apply the method of QSC (Quantum Simulated Cooling) to the Hamiltonian ana- 
logue of Grover's Algorithm [24]. Although the arguments given in this section are 
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heuristic, the claims made for the general scheme are rigorously shown in the sup- 
plemental. Say that we are given a function / : {0, 1}" — > {0, 1}, and we wish to 
find y g {0, 1}" such that f{y) = 0. Analogously, say that we are able to simulate 
a Hamiltonian Hs on n qubits, which has only two distinct eigenspaces. These are 
labeled Pq and Pi , with energies = and wi , and correspond to the logical qubit 
states \x) such that f{x) = (1), respectively. For notational simplicity, we let the 
symbol for an eigenspace also represent its projector. We then have 

Hs - c^iPi (2.1) 

Up to a constant factor, the simulation of Hs for a fixed time is equivalent in cost 
to evaluation of f(x) [26], so solving for y is equivalent to finding a state in the zero 
energy manifold of (2.1) [27, 26]. 

To take the system from an initial state \F) to the ground state manifold, we 
concurrently evolve a single qubit bath, with Hamiltonian Hb = i^B^s ® lt)(t|j and 
define the non- interacting Hamiltonian H ~ Hs + Hb- We prepare the system and 
bath in the state \F and introduce a coupling between them, denoted by 

V = no\F){F\ ® (2.2) 

To prevent accidental symmetries leading to frustration, \F) is a randomly generated 
quantum state [28]. 

To illustrate the evolution of our state with this interaction, we decompose 
into its spectral components: 

|F) =a;o|0)+xi|l) (2.3) 

where Xj \j) ~ Pj \F) and Xj is real and non- negative. As the space S = 
Span{|0 \.) , |0,t) , jlji) , jilt)} is invariant under H + V, if we prepare the system in 
the state \F), we may reduce our analysis to this subspace. Written explicitly, within 
S the full Hamiltonian takes the form 
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Suppose that we set fio ^ wi, and lob ~ ^i- Since |0 f) and |1 \.) are nearly 
degenerate, we expect that superpositions of these states form eigenstates oi H + V. 
As V has no diagonal terms in the above basis, we see that even orders of flo induce 
level shifts on these states, while odd orders couple between them. With knowledge 
of xq and xi, we could compute [29] the even-order level shifts induced by V on the 
energies of |0 t) and |1 t), and adjust lub accordingly so that they are degenerate. 
To leading order, determining lob is equivalent to solving ujb H — ^j-^ = oji -\ — ^j-^- 

Notice that for \xo\ <^ \xi\ ^ 1, the level shift of state |1 ^) is approximately ^J' , 
which is non- negligible compared to the coupling floXoXi. As seen in Figure 1, if this 
shift is not accounted for the scheme's success rate becomes exponentially small with 
increasing n. 

Since the interaction strength JIq is perturbative, in the limit \xo\ <^ 1 we conclude 
that for the manifold of states with energy near lub, the Hamiltonian is effectively 

i/e//=C.B(|Ot)(Otl + IU)(UI) 

+r!oXoa;i(|Ot)(U|+h.c.) 
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Fig. 2.1. This figure illustrates the sensitivity of the Grover cooling scheme to errors in bath 
detuning, as a function of n. As seen in the inset for n = 7, the probability of the cooling scheme 
to produce the ground state can be computed as a function of bath detuning, {ljb — i^i)/<^i- The 
vertical dashed line represents the first order correction to the level shift between states |1 J,) and 
|Ot}j which is Qq/uji for \xi\ 1. The main figure plots the half-width of this fidelity function, 
for Qq/uji = 0.05. The dashed line represents a linear fit, suggesting that the half-width scales 
approximately as \/2~"aJi/f2o- Naively setting uib = ^1 requires that, for a fixed success probability, 
flo/uji must be scaled as \/2~", thereby negating the quadratic speedup observed with the optimal 
detuning, 

while the states |0 I) and |1 t) are unchanged hy H -\-V. We thus observe coherent 
oscillations between |0 f) and |1 |), at a rate = ^ox^Xj. If we prepared the system 
and bath in the state i), we could evolve for a time t = The resulting output 
would be 

^-^r{H+V) (^^ |o 1^ ^ |i ^ |o ^ |o 

SO that the system is in a groundstate of Hs ■ 

We can relate the actual cost of the algorithm to the simulation time by noting 
that a single implementation of exp(— itiJs) is equivalent in cost to 0(1) evaluations of 
the function / [26]. Using a stroboscopic expansion [30], simulation of Hs+Hb + V for 
total time T may then be done on a standard quantum computer at a cost approaching 
0(||if||T) [31]. Since the time scale necessary to map between |1 \) and |0 t) is set by 
the Rabi rate Q. = Q.qXqXi, one sees that for fixed wi, the total cost of the algorithm 
scales linearly with |a;o|^^- The scaling of simulation time as |xo|^^ will also apply to 
generalizations of QSC to more complicated Hamiltonians. 

Say that iV = 2" is the dimension of the ti— qubit Hilbert space, and iVo = 
dim Po- If we select \F) from a random sample, so that on average |a;o|^ ~ Nq/N, the 
running time of the algorithm will scale as -^Z N/Nq , reflecting the quadratic speedup 
over classical computation observed in Grover's Algorithm. On a standard quantum 
computer, such a sampling can be achieved through e-approximate unitary 2-designs, 
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Fig. 2.2. The spectrum of Hg + Hb during the cooling of energy level j. The left and right 
columns represent the spectrum of the composite system when the bath spin is in its ground and 
excited state, respectively. Each band (dashed line) represents a degenerate subspace of Hprop, 
whose degeneracy is lifted by Hinp^t ■ The bath gap ujg is adjusted so that the ground state band on 
the right is nearly degenerate with band j on the left. The inset (center column) shows the effect of 
the interaction V on the two unperturbed eigenstates in Si, \j 4-} and |0 t)- The level shift in these 
states is due to even order corrections in V, and the value of u)g is adjusted so that these shifted 
states are exactly degenerate. The odd order corrections in V then allow for coherent oscillations 
between these states at frequency 2f2. After interacting for time 7r/(2f!), the bath spin is pumped 
back to its ground state and the process is repeated for band j — 1. 



which may be implemented at a cost of 0(nlog(l/e)) [32]. Note that being able to set 
T correctly, as well as correcting for the level shift induced by V, requires knowledge 
of the value of xq and xi. This issue is also relevant to the more general problem, 
and in the case where the decomposition of \F) is unknown, we present a modified 
scheme below that succeeds probabilistically in the same time. 

3. Cooling a Quantum Circuit. Here we show how QSC can be used to pro- 
duce the outcome of a chain of 2-qubit unitary operations, U = UlUl-i---Ui, at a 
cost scaling as 0(poly(L)). Since 1 and 2-qubit unitaries are sufficient to implement 
any efficient quantum computation [33, 34], any problem efficiently solved through 
standard quantum computation can also be solved using QSC with at most a polyno- 
mial overhead. The idea behind our result draws from the work of [35], which shows 
that adiabatic quantum computation is equivalent to standard quantum computation. 

Suppose there exists a Hamiltonian Hs whose unique groundstate, after tracing 
out any ancilla qubits, can be made arbitrarily close to U |0"). Preparing the ground- 
state of Hs would correspond to producing the outcome of the computation. One 
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Hg satisfying this requirement is a variant of Kitaev's eloek Haniiltonian [36]. As in 
[37], to deseribc Hg we consider a particle living on a ID lattice with L + 1 sites, 
whose internal state is described by n qubits. For a given site the particle has fixed 
onsite energy w, but may also tunnel to neighboring site I + 1 through a coupling term 
— ^■Ui acting on the internal states. The Hamiltonian describing the particle is then 

L ^ L 

Hprop/W = \si){si\ - 2 I] {Ui\si){si^i\+h.c.) (3.1) 

(=0 (=1 

where \si) corresponds to the particle being in the Ith site. This Hamiltonian is 
analogous to that of a particle freely propagating through space, and its eigenstates 
are all of the form [38]: 

L 

Co la;) |so) + ^ Q {UiUi^i...Ui) \x) \si) (3.2) 
1=1 

where \x) is the internal state of the particle at site I — 0. 

Since we want |0") at the start of the computation, we add a perturbation to 
Hprop of the form 

n 

in 

(g) \so){so\ (3.3) 

?7l — 1 

where Ai <C and |l)(l|m acts only on qubit m of the particle's internal state. This 
lifts the degeneracy between eigenstates with different initial conditions, and allows 
us to reduce our analysis to the invariant subspace of states with \x) = |0") in (3.2). 
We label this subspace Si and its complement 5*2. 

The ground state of Hprop in Si corresponds to c/ = -yj^ for all I. The site 
L component of the ground state is "^^|0") I^l), so by preparing this state and 
measuring the particle in L we would obtain the outcome of the computation. If we 
instead concatenate 0{L/e) identity operations to the definition of Hprop, then after 
tracing out the particle's position its internal state has an 0(e) trace-norm distance 
from J7|0")(0"|J7^ [35]. 

The eigenstates of Hg = Hprop + Hinput in Si are non-degenerate and have energy 

cjj = cj • ^1 — cos (^-g^^^ , where < j < L. There are two other important energy 

scales associated with Hs- The first is the maximal energy, ||^^s|| = w • 0{L) [35], 
where 1 1 • 1 1 is the spectral norm. The other relevant energy is the spectral resolution 
of ^i: 

A = min{|£; - ujA : E e Spec(i7s), E ^ ljA (3.4) 

As seen below, the ratio ||i?5||/A, where || • || is the operator norm, will determine 
the overall cost of the simulation. Further, for a given error tolerance the energy scale 
A provides an upper bound to the system-bath coupling V and allowable error terms 
in the sumulation. 

To bound A, we note that the eigenspaces of Hprop are degenerate bands corre- 
sponding to the spectrum {i^j}j"=0' which are indexed by the state \x) in (3.2). Hinput 
is diagonal within each ujj band, with diagonal entries AiN{x)hj , where N{x) is the 
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number of I's in the binary expansion of \x) . By diagonalizing Hp^op [39], one may 
verify that 1/hj ~ 0{L^). Sinec Hinput vanishes only when \x) ~ |0"), as long as first 
order degenerate perturbation theory is valid we see that the spectra of Hs in Si and 
^2 are distinct. This is true when ||i?i„put|| = nl^i is bounded by the minimum level 
spacing of Hprop, which scales as oj ■ L~^. We therefore assume that Ai cx w • ri^^i^^, 
so the gap satisfies A"-'^ < (Aihj)^^ = uj~^ ■ 0{nL^). 

To prepare the groundstate of Hs, we simulate a single qubit bath coupled to 
the system, with energy splitting lob- The projectors into Si and 5'2 act trivially 
on this qubit. The fiducial state is \F I) = |0") |so) \i} & Si. We also introduce the 
system-bath interaction V = ftoln <E) |so)(so| ^ (Jx- Like Hs, V is block diagonal in 
Si and S2 and satisfies 

SiVSi=^lQ\F){F\®ax (3.5) 

Letting \ j) denote the j-th excited state in Si, we write the spectral decomposition 
of 

\F)=Y.^,\j) (3.6) 
3 

Although the states |j) are dependent on U, the coefficients Xj are only dependent 
on L and may be calculated explicitly [39] . 

The algorithm proceeds as follows. We start with the bath energy ujb near lol, to 
address the transition \L \) — > |0 t) (see Figure 2). In order for this to be favorable, 
we must calculate the even-order level shifts induced by V on states |0 f) and \L I), 
and adjust ujb so that they are degenerate. Such a calculation is equivalent to finding 
a root of a degree L polynomial to accuracy l/poly(L), and is possible since one may 
explicitly calculate the coefficients Xj in (3.6). We evolve the system and bath under 
the Hamiltonian iJ -|- y for a time tl = where VI = ilol^^oa^LKl + O(51o/^)) is 
the coupling rate between |0 f) and \L \). This maps the component of \F \) in the 
state |L I) to the state [Of), while leaving the lower energy space Span { | j ^) } ^""0^ 
unchanged. We then measure the bath in the logical basis. A measurement of jf) 
implies that the desired transition has occurred, so we may terminate the algorithm. 
A measurement of ||) implies that we have projected to Span{|j ■ We therefore 

decrement ujb to a value near ujl-i, account for the level shifts of |0 f) and |(L — 1) J,), 
evolve for time ri_i, and remeasure the bath. Repeating this process for at most L 
evolutions, we reach the ground state with high probability. 

To prove this claim, we use the language of trace-preserving, completely positive 
(TCP) maps [40, 41, 42]. The map associated with measurement and conditional 
evolution at energy near ujj is 

E,{p) = u,\m\p\m\u] 
+ it)(tiPit)(ti 

where p is a density matrix for the system and bath, Uj is the unitary evolution 
under H + V ior ujb ~ and time r,- « Wh^^^^ ^^'^ = Is ^ |i)(i|B- 

We define the TCP map describing the complete algorithm as E = i?i o £'2... o El- 
In the supplemental section, we build upon results from [38, 43] to show that for 



E(p) = Aopo + R 



(3.8) 
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where Ao > and pq is a density matrix made up of states |0 l) and |0 f)- R represents 
an error term satisfying Tr[i?] = 0(L'^/^fio/A) = ||^||. For the algorithm to succeed 
with error rate e = Tr[(l - |0)(0|)E(p)] = Tr[7?,] - Tr[|0)(0|A], we must scale as 
0(eAi-5/2). 

3.1. Timing and Cost. We now discuss the relationship between the simulation 
time T and the actual cost of implementing the algorithm. We assume that the cost 
of simulating n qubits under a Hamiltonian H for time T scales as poly ||T). 
This is true for all H that may be written as a sum of poly(n) interactions, each 
involving a fixed number of qubits [2] . For our clock Hamiltonian, we may implement 
Hs + Hb + V using at most 5-qubit interactions [35], and expect that a 2-local QSC 
scheme should also be possible using methods in [44, 45]. 

To leading order in (ilo/A) the time required for cooling step j is Tj — boa; I ' 
so the total time of the algorithm may be computed as 



By diagonalizing Hprop within the space Si, a simple calculation shows that the coef- 



ficients Xj in (3.6) satisfy \xj\-^ = O cos (^^^f^ j j . Since flo = 0{eAL-^/^), 

the total simulation time T scales as A^^ • O (^L^\/Llog{L)/ej . Finally, as A = 
0{n^^ L^^) and ||iJ|| = 0{L), the actual cost of the algorithm is then set by = 
O (^nL^'^^/L\og{L) /ej . QSC is therefore equivalent to standard quantum computa- 
tion: any problem on n qubits that may be solved in poly(n) time using a standard 
quantum computer may also be solved in poly(n) time through QSC. 

4. Extension. Although the scheme above is in principle as powerful as circuit- 
based quantum computation, its applicability towards cooling other Hamiltonian sys- 
tems is limited by the assumptions placed on Hs and V. To avoid unwanted tran- 
sitions, it requires an accessible non-degenerate eigenspace Si with gap A, as well 
an interaction of the form (2.2) described by a known decomposition (3.6). Here we 
propose an extension to this scheme that overcomes these constraints. The caveats 
of this approach are that it is probabilistic and requires {k -\- l)-local interactions for 
fc- local Hs- 

We briefly summarize the scheme before describing it in detail. We assume that 
the spectrum of Hs has a non-degenerate ground state |0), as well as a manifold Pi of 
energies near wi that is well separated from the rest of the spectrum. We also assume 
the ability to simulate evolution under to an operator Ts that couples between this 
manifold and |0). We define Hamiltonians H and V based on Hs and Ts such that 
coherent oscillations occur between |1C) and \0B), where |1) £ Pi and \B),\C) are 
states of a qutrit bath. We start with a fiducial state \FC) and evolve under H + V 
for a sufficiently long time, after which we check if the bath is in state \B). If that is 
the case, we verify that the system has transitioned to the ground state by mapping 
\B) — >■ \L) and evolving under a verification Hamiltonian (4.8). A bath measurement 
of \R) heralds the success of the scheme, and otherwise we repeat the process. 

We now specify the scheme's requisite assumptions. First, the Hamiltonian Hs 
is of the form 




(3.9) 



Hs = PiHsPi + P2HSP2 



(4.1) 
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where Pi and P2 arc projectors into eigenspaces of the same name. The eigenspace 
Pi represents a narrow band of energies within (uji — Suji, uji + Suji). We assmnc that 
Hs has a non-degenerate ground state |0) with energy wq = 0, and let P2 represent 
the space orthogonal to both Pi and |0). The limiting energy scale for this scheme is 
A = min{a;i,i?, |i? — wi| : _E G Spec(i/5|p2)}. A^^ will set an upper bound for the 
time scale of the simulation. 

Along with Hs, we assume access to a Hermitian operator Ts that couples |0) to 
the space Pi: 

PiTs\0)=n\l) (4.2) 

where is real and positive by choice of phase convention. Using this coupling, 
over a time scale l/fi the simulation will cause coherent oscillations of the form 
|1C) -O- \0B), where |C) and \B) are orthogonal bath states. We must further assume 
that the energy spread of Pi is small compared to the coupling: 6u!i <^ fl. This will 
allow us to treat |1) as an eigenstate of Hs in our analysis of the simulated evolution. 
We must also have that the energy shift induced by T5 on the ground state is small 
compared to the coupling: | (0| T5 |0) P/wi ^ ^, that Tg is perturbative compared to 
the gap: 51 < IIT5II ^ A, and that T5 does not couple strongly between |0) ,Pi and 
P2: ||(|0)(0|+Pi)T5P2||« VHA. 

The original scheme succeeds by introducing a 2-level bath and adjusting its en- 
ergy ujb so that |1 4,) and |0 f) are nearly degenerate. Within degenerate perturbation 
theory, the interaction V then produces a splitting of 251 between approximate eigen- 
statcs (1 1 i) i 1 t) ) ! thereby causing coherent oscillations of the form 1 1 1) o 1 t) • 
If not accounted for, the (even-order) level shifts induced by V on the state |1 1) can 
be much larger than this splitting, meaning the eigenstates look more like |1 ],) and 
|0 t), so that the desired oscillation does not occur. Although we may explicitly ac- 
count for the level shifts by adjusting ljb, this requires knowledge of the coefficients 
Xj in (3.6). Instead, one may tailor the unperturbed Hamiltonian H so that unwanted 
level shifts cancel out: 

H = Hs^ (|C)(C| + \R){R\ - \L){L\) 
+ ujils {\R){R\ + \L){L\) 

The bath Hilbert space now has dimension 3, with basis vectors \C) , \R) and \L). 
We use the operator T5 to create the system-bath interaction V: 

V = Ts®{\C){B\ + \B){C\) 

1 (4.4) 

\B)^—i\L) + \R)) ^ > 

To see how H + V causes the desired oscillations, we compute the level shift operator 
I](z) for the manifold P of eigenstates of H with energy at most A/4 away from uji 
[29]. Written as a projector, 

P = |0)(0| <E) i\L){L\ + \R){R\) +Pi<E) \C){C\ (4.5) 

E(z) characterizes the time evolution of the states in P under Hamiltonian H + 
V. As seen in [35] and further developed in the supplemental, the eigenstates and 
eigenvectors of S(z) for z close to uji are good approximations of eigenstates and 
eigenvectors oi H + V. Written explicitly, we have 

E(z) = PHP + PVP + PVG{z)VP , , 

(4.6) 

-^PVG{z)VG{z)VP + ... ^ ' 
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where G{z) = ^qJ^hq ^^'^ Q complement of P. We observe that for z = Ui, 



+ o.r'|ocxoci + ^\ » icxci 

so that {B\ G{uJi) \B) = 0. Thus (4.6) truncates at second order for z = oji, so that 
the only contribution from V is PVP = fl {\1C){0B\ + |0B)(1C|) and the ground state 
shift lMkli>l!|0B)(0S|. As seen in the supplemental section, H^jf = is a good 

description of the evolution of states in P, and by our bounds on ' J^' and Soji, 
it causes the desired oscillations at a rate 2n. 

Say that our fiducial state \F C) has overlap |/i| = Since coherent 

oscillations between |1C) and \0 B) occur at frequency 251, as long as we evolve 
for times sampled randomly over a range ^, with probability 0(|/ip) we expect to 
observe a transition of the form |1 C) — )■ \0 B), heralded by a measurment of the bath. 

The average simulation time of the algorithm then scales as O (j/Tpn) • If o^i^ is not 
given an explicit value of f2, as in [46] one may implement the scheme with evolution 
times sampled randomly from [r, 2t] , and iteratively increase r — >■ 2t after ~ 1/ 1 /i | 
failed attempts. Since the sampling time r grows exponentially, the total evolution 

time before success still scales as O ■ 

Unfortunately, a bath measurement of \B) does not imply the system is in its 
ground state, as other resonant transitions could also occur. To account for this, after 
measuring \B) we map the bath state \B) — > \L) and evolve under 

H + nols®{\L){R\ + \R){L\) (4.8) 

for time r„ = . The inversion symmetry of H implies that only the system ground 
state |0) has degeneracy between its \L) and \R) bath states. All other eigenstates 
lip L) have at least A less energy than any \R) eigenstate. Hence by energy conser- 
vation only |0) exhibits coherent oscillations from \L) to \R), so measurement of the 
bath in \R) heralds success of the scheme. 

Notice that, as long as 5coi <^ A, all of the constraints on Ts required for this 
scheme are satisfied by Ts = ilo\F){F\. Since \F){F\ is rank 1, the assumption that the 
groundspace Pq is non-degenerate is no longer necessary, as Ts vanishes on every state 
in Po orthogonal to /o |0) = Pq \F) (where \fo\'^ = \ {F\ Po\F) \). In this case we would 
have = nofoPi \F), so [/ip = {F\P-l \F) and n = |f^o/o/i|- If we use unitary 
2-designs to randomly generate \F), we require that there is a fixed probability in n 
that n > flo^/do ■ di/2"+i and l/ip > di/2"+i, where di is the rank of Pi. Thus if 
we are given no information about Hs other than Jwi, wi, and A, by using 2-designs 
and the probabilistic scheme we may obtain the ground state of Hs with an average 
simulation time scaling as 0{ ^^jl^x/^ ^ ), reflecting the quadratic speedup observed in 
Section 2. 



5. Concluding Remarks. There are several known alternatives to standard, 
logic-based quantum computing [22, 47, 48, 49, 50, 51, 52]. The advantages of our 
scheme are that it requires the simulation of only time-independent Hamiltonians, as 
well as measurements of a single qubit (or qutrit) bath. Using the gadget construction 
[44, 45] we expect that it may be efficiently implemented with only 2-local interactions. 
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Our work suggests several new avenues for investigation. One question is whether 
the techniques used in QSC may be apphed to prepare other interesting states, such 
as mixed state ensembles [10, 53, 54, 55] or ground states of frustration free Hamil- 
tonians [51]. Using techniques in [56, 57, 58], one may attempt to show whether this 
scheme is robust against time dependent error terms in the simulation, or nonuni- 
tary evolution described by weak interactions with an environment. Finally, we note 
that the application of QSC to the clock Hamiltonian in Section 3 used only a 2- 
body system-bath interaction and a product fiducial state. Similarly to topological 
quantum computing, this prompts the question of when local interactions suffice to 
produce the ground state of a Hamiltonian, and fundamentally, what the relationship 
is between a Hamiltonian's computational complexity and the potential to cool it 
using such interactions. 
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Supplemental Sections- 

Bclow we make rigorous the claims stated in the main body of the work. We 
start by developing some preliminary mathematical tools, then go on to give suffi- 
cient conditions for the success of the deterministic and probabilistic (extended) QSC 
schemes. We conclude by analysing the modified Kitaev clock Hamiltonian, thereby 
showing that both forms of QSC are polynomially equivalent to standard quantum 
computation. 

6. Mathematical Tools. The following theorems ensure that the effective Hamil- 
tonian, Heff, that we derive accurately describes the dynamics of the simulator. The 
first result is a slight modification of Theorem 3 in [38]. The theorem is concerned 
with a Hamiltonian H and a perturbation V to the Hamiltonian. Within the sub- 
space of interest, the theorem gives a one-to-one correspondence between the spectra 
of Hpff and H = H + V , and bounds their difference. Before stating the theorem we 
require a few definitions. 

Let M = P(BQ describe a finite dimensional Hilbert space on which H acts, where 
P is spanned by the eigenvectors of H whose eigenvalues are in (A_, A-|-) . Likewise 
define P, Q with respect to H , using the same bounds. For simplicity we let the 
symbol for a subspace also represent its projector. We will assume that H has gap A, 
i.e. that the energies of H in P are at least A away from those in Q. We are interested 
in the dynamics under H\p, which we approximate with H^ff- The approximation is 
derived from a series expansion of the self-energy operator [29] : 

Sp(z) = PHP + PVP + PVQ{z ~ Q{H + V)Q)-^QVP 

= PHP + PVP + PVQG{z)QVP + PVQG{z)QVQG{z)QV P + ... (6.1) 

In this notation G{z) = (2 — H)^"^ is the Green's function for the unperturbed Hamil- 
tonian. 

In all cases below, || • || represents the operator 2-norm, 

sup \\X\v)\\ 

{v\v) = l 

where X is a (bounded) linear map from H to a finite Hilbert space H'. and || \ v) \ \ 
is the norm induced by the inner product of H'. This is a consistent norm, satisfying 
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< \\A\\ ■ \\B\\ [43]. Finally, wc mention a slight abuse of notation: If |f) is a 
vector in a Hilbert space and A an operator acting on that space, then for expressions 
of the form 

\v)+0{r) 
A + 0{r) 

0{r) represents a vector (operator) with norm scaling as 0(r). With this notation 
in hand, we can state the first theorem. Note that unless otherwise mentioned, the 
proofs for the following results are at the end of the section. 

Theorem 6.1 ([38]). Assume that H has no eigenvalues in [A_ - A/2, A_ + A/2] 
and [A_|_ — A/2, A+ + A/2], and that \\V\\ < A/2. Assume that there exists an operator 
Heff on P whose spectrum is contained in [c,d], and that for some ^ > Q, we have 
that 

[c-7,d + 7] c (A_,A+) 

\\T.p{z) - HeffW < 7 

for all z ^ [c — 7, d + 7] . Then if Xj (Aj) is the jth largest eigenvalue of H^ff {H\ p), 

\ <7 

The result of [38] is for the case A_ = —00. Since the proof of Theorem 6.1 only 
requires a straightforward modification of the original, we do not show it here. The 
following result is used in the proof of Theorem 6.1, as well as in some of the claims 
below. 

Lemma 6.2 ([38]). Let H, H he two Hamiltonians with ordered eigenvalues 
Ml < M2 ^ •■• o,i^d (Ti < (72 < ... • Then for all j , 

\fij-a,\<\\H-H\\ 

Theorem 6.1 gives bounds for the error in the approximate eigenvalues of Hf,ff, 
but in order to sufficiently describe the dynamics we also need a correspondence be- 
tween the eigenvectors oi H and Heff- To that end, we give a result derived Theorem 
3.6, Chapter V, of Stewart and Sun's Matrix Perturbation Theory [43]. It is effectively 
a statement of the conservation of energy, and will ensure that transitions which do 
not preserve energy are suppressed. 

Theorem 6.3 ([43]). Let H, H be Hermitian operators. Suppose H is resolved 
by P and Q: H = PHP + QHQ, and 

Spec{H\p) C [A_ + A/2, A+ - A/2] ^ 



Spcc{H\q) C (-00, A_ - A/2] U [A+ + A/2, 00) 



12 



Algorithmic Cooling of a Quantum Simulator 



Let P be the span of the eigenvectors of H with eigenvalues contained in (A_,A_|_). 
Then for any \v) G P, \v) G P, 



(v\P\v) > 1 - 



2\\H -H\ 



A / 

The proof of Theorem 6.3 requires the foUowmg lemma. 

Lemma 6.4. Let A, B he Hermitian operators on C^jC", respectively. Suppose 
that \\A\\ < a and that B is invertihle, with < ^jqrg, for /3 > 0. Then for any 

linear operator X : C" — > C"\ 



The proof of Lemma 7.2 in the following section requires the following result. It 
will be used to bound the effect of error terms during time evolution: 

Lemma 6.5. Let A, B he Hermitian operators on some space Si® 82- Suppose 
that A = SiASi + S2AS2, and that 

\\SiA~^Si\\ < l/Gi 

\\S2A-'S2\\ < I/G2 
foil = \\SiBSi\\< Gi/2 
b22 = ||'S'2-B5'2||< G2/2 
612 = \\SiBS2\\< min(Gi,G2)/2 



Then (A — B) is invertihle, and 

\\SM-B)-^S2\\< 



612 



(Gi -6ii)(G2-622)-6L 



\\Si{A-Br^Si\\< ^ ^ , (l + hi2\\Si{A-B)-^S2\\) 
Gi — Oil 

\\S2{A - B)-^S2\\< ^ \ (1 + bi2\\SM - B)-^S2\\) 

Lr2 — O22 

The goal of these theorems is to bound the error in the unitary evolutions de- 
signed to map \i 4,) — >■ |0 t) or |1 G) — > |0-B). This will be done by showing that both 
the eigenstates and eigenvalues corresponding to Heff are close to true eigenstates 
and eigenvalues of H . Theorem 6.1 is used to characterize the spectrum of Heff - The 
following corrolary states that if an eigenspace P' C P of H^ff is 'well resolved' from 
its complement, then the corresponding eigenspace of H is well approximated by P'. 
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Corollary 6.6. Given the assumptions of Theorem 6.1, let P' <Z P be an 
eigenspace of H^f / , and let P — P' be its local complement. Define: 



V 



= max{|a; - y| : a:,y G Spec(i?e// |p')} 



77 = min{|a- -y\:x<^ Spec{Heff\p'),y S Spec(iJe//|p-P')} 

Define P' as the eigenspace of H obtained from the eigenvalue correspondence in 
Theorem 6.1. If rj > 7, then for each eigenstate \vi) € P' of H, 




P' \v,) > 1 - 



Corollary 6.7. Let H be a Hamiltonian resolved by spaces P and Q: H = 
PHP + QHQ. Assume that P ^ Pi® P2® Pl, and that there exist < Ai+ < 
A2- < ... < Xl+ such each Pk corresponds to the eigenspace of H with energies in 
[Afc_ + A/2, Afc+ — A/2]. Further assume that the eigenvalues of H in Q are at least 
A > away from those in P. Finally, say that Pk is the eigenspace of H corresponding 
to eigenvalues within {Xk-,Xk+), for I < k < L 

Then for each \v) G P, \v) G P = Pi © ■•■ © Pl, 




We now begin the proofs of the above results, neglecting Theorem 6.1 and Lemma 
6.2 as they may be derived (with minimal modification) from results in [38]. 

Proof of Lemma 6.4-' Since )] • )) is a consistent norm, we have that 

\\AX\\ < a\\X\\ 



\X\\ = \\XBB-^\\ < 



a + l3 



\\XB\\ > {a + p)\\X\\ 
By the triangle inequality we conclude that 

\\AX - XB\\ > \\XB\\ - \\AX\\ > I3\\X\ 



□ 
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Proof of Lemma 6. 5: 

First, we show that [A — B) is invertible. It suffices to show that < 1, as 

then (1 — BA^^) — ((A — B)A~^^ is invertible. Let a norniahzcd \vi) G Si be given. 
Then A^^ \vi) = c\v[) for some (normahzcd) \v[) G Si and |c| < l/G^. Using the 
above equation, one gets 



\BA-^ \vi) IP = {vi\A-^BBA-^ \v. 
1 



<^{v[\B{S, + S2)B\v[) 



- G1 
< 1/2 

For any normahzed state \v), one may write \v) ~ a \vi) + b \v2), with \vi) G Si and 
|ap + |6p = 1, and by the triangle inequality 

\\BA-^ \v) II < |a| • \\BA-' \v,) \\ + \b\ ■ \\BA-^ \v2) \\ 
<±{\a\ + \b\)<l 



where the last inequality follows from writing |a| = ^/s, \b\ = -^/l — s for some s G [0, 1]. 
This shows that ||_BA~^|| < 1, which implies that (1 — BA~^) has only positive 
eigenvalues and is therefore invertible. 

Since we have shown that (A — B) is invertible, one may easily check that 

{A- B)-^ = A-^ + A-^B{A- B)-^ (6.2) 

Decomposing into Si and 5^2 components, one gets 

Si{A - B)-^Si = SiA-^Si + SiA-^Si {SiBSi{A - B)-^Si + SiBS2{A - B)-^Si) 
By the triangle inequality and the relation HCi^ll < ||C|| • \ \D\\, we have 

\\Si{A-B)-'Si\\ <^ + ^\\Si{A-B)-'Si\\ + ^\\SiiA-B)-^S2\\ =^ 
\\Si{A~B)-^Si\\ <—±—{l + bi2\\SM-B)-^S2\\) (6.3) 

Grl — Oil 

where we also used the fact that we may interchange the projectors Si, S2 when taking 
the operator norm. Interchanging the numbers 1 and 2, we obtain an equivalent result 
for ||S2(A — £?)^^S'2||. This proves the last two statements in Lemma 6.5. 
Likewise, using (6.2) one may show that 

||S'i(A - B)-'^S2\\ = \\SiA-^Si {SiBSi{A - B)-^S2 + SiBS2{A - B)-^S2) \\ 

< ^-21\\SM-B)-'S2\\ + ^-^\\S2{A-B)-^S2\\ =^ 

(_rl Grl 

\\SM - Br^S2\\ < „ \\S2{A - B)-'S2\\ (6.4) 

Oi — Oil 

Substituting the result of (6.3) (with 1 and 2 interchanged) into the right hand side 
produces the first statement of Lemma 6.5. 
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□ 



Proof of Theorem 6. 3: 
Define the numbers 



2 

Note that A > A/2 since (A+ - A/2) > (A_ + A/2). Since H is Hermitian, there 
exists a unitary operator U ~ {Xp,Xq) that diagonaUzes H, where the columns of 
Xp and Xq form an orthonormal basis for P and Q, respectively. We see that 

H — \\ = XpLpXp + XqLqXq 

where Lp and Lq are diagonal matrices, with eigenvalues in [—(A — A/2), A — A/2] 
and (— oo, —A — A/2] U [A + A/2, oo), respectively. Analogously, we may define the 
decomposition oi H — X: 

H-\1 = XpLpXI,+XqLq^ 

where the eigenvalues oi Lp and Lq are in [—A, A] and (— oo, —A] U [A, co). 
Consider the operator 



VE = XliH-H\XQ 



From the identities above it follows that 



Ve = xliH - X1)Xq - X\{H - X1)Xq 



p\ J p\ 

= LpX^-Xq - X^pXqLc, 



p W p 

Noting that ||ip|| < A and ||Lq^|| < ^^^y^ , we may use Lemma 6.4 with A = Lp, 
B = Lq, a = A, P = A/2, and X = X^Xq to conclude 



_ 2||xt (^H-h)Xq\\ 

~ A 

2\\H~ H\\ 

Where the last line follows from < ||A|| • 

Let \v) £ P he given with {v\ v) — 1. Since the columns of form an orthonor- 
mal basis for P, there exists a vector \x) such that \v) ~ Xp\x) and {x\x) = 1. From 
the previous line we conclude that 

{v\Q\i)^{v\XQXl\i) 

= {x\xIXqXIXp\x) 
t 



< WX^Xs 



pi 



2\\H -H\ 
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where in the last hne wc used the fact that the operator 2-norni is unchanged when 
taking the Hermitian adjoint. The first statement follows by noting that Q = 1 — P. 
To prove the second statement, we can make an identical argument using Ve = 
- H)Xq. 

□ 

Proof of Corollary 6. 6: 

By assumption, we have that H = PHP + QHQ, where P projects into the 
eigenspace of H with energy in (A_ + A/2,X+ — A/2). This defines the energy gap 
between H\p and H\q. For H = H + V, we define P as the eigenspace of H with 
energy in (A_, A+). P' C P is an eigenspace of Hf,ff and P' C P is the associated 
eigenspace of H, obtained from the eigenvalue correspondence discussed in Theorem 
6.1. 

Let \vi) £ P' be an eigenstate of H, with corresponding energy Ei. We may write 

= N, \u,) + Q 

where P \vi) — Ni \ui). Since \vi) S P, we may use Theorem 6.3 to conclude that 



A 

Now decompose \u.i) into components parallel and orthogonal to P': 

\u,) = P' \u,) + {P- P') \ih) 

We wish to bound {ui \ P' \ui) from below. To do this wc first note that, as in the 
proofs of Lemmas 5 and 6 of [38], Y,p{Ei) \ui) — Ei \ui). From this we conclude that 

(I]p(^,) - H^ff) \u.,) + {H,ff - Ea)P' \u.d = -{Heff - Ea){P - P') \u^) 

By assumption, \\{Y,p{Ei) — iJe//)|| < 7. Since Ei is an eigenvalue corresponding to 
P', the eigenvalues of operator {Heff — Eil)P' have magnitude at most i/, and since 
\Ei — Ei\ < 7, \\{Heff — Eil)P'\\ < J + v. By the triangle inequality, the norm of 
the left hand side is less than 27 + f. Likewise, the eigenvalues of [H^ff — Eil)\p^p> 
have magnitude at least rj so the eigenvalues of {Heff — Ei\)\p-pi are greater than 
77 — 7 > 0. Thus the norm of the right hand side is greater than (77 — 7)||(P — P') \uij \\. 
Therefore 

ll(P-P')k)ll<^ 
77-7 



27- 



?/ - 7 



{u,\ P' \u,) ^ 1 - {u,\ (P - P') \u,) > 1 - 

The result follows by noting that (vi\ P' \vi) = Nf {ui\ P' \ui). 

□ 

Proof of Corollary 6.7: 

We may assume without loss of generality that each domain + A/2, Afc+ — A/2] 
is at least 2A away from its neighboring domains. If this were not the case for a 
neighboring domain, the gap A from Spec(i7|Q) would imply that there exist no 
QHQ eigenvalues between the two domains, so they may be merged. 
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Let a normalized vector \v) e P be given. We may decompose \v) into its com- 
ponents within each P^,: 

where J2k ~ ^- a-Pply Theorem 6.3 to any given P^ to get 

l-\rk\^ = {vk\P\vk) > {vk\Pk\vk) 
where \rk\ < HiliL_£ii^ go for each fc, 

where 

Substituting into the definition of \v), we get 

\v) ^Y.'"' {p\^''^^+'^'^K)] 

k 

^p(^Ck\vk)^ +reff\v^) 

where r^ff l^^} — '^k'-^fk \''^k) ^ vector in Q. By the triangle and Cauchy- 
Schwartz inequalities we conclude that 

k 

'2\\H -H\\\ 



k 



k \ 

' %\ii-R\ 



A / 

The first statement of the corollary follows by noting that (v\P\v) = 1 — |fe//p. The 
second statement follows by making the same argument and, except with the initial 
assumption, adding or removing ^ to each projector and vector. 

□ 

7. Deterministic QSC Analysis. This section describes the deterministic Quan- 
tum Simulated Cooling scheme, and gives sufficient conditions for its success up to 
an infidelity 0(e). 

Theorem 7.1. Let Hs he a Hamiltonian acting on n qubits, with eigenspaces Si 
and S2- Assume that 

L 
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where Q — ojq < uji < uj2--- < o,nd 

A min{|£' - lo^\ : E <E Span(iJs), E 7^ cjj > 

Ui,E 

Finally, assume that there exists operators Ts, 5j, such that 

SiTsSi = r!o|G)(G| 



with 



L 

\G)^Y.'xA^) 



\xj\ > 0, and for all j , 



r = no/ A < 1/8 

< r • riolccoXjI/A 



|S'i(5jS'i| 



A 

\Si{Ts + 6,)S2\\^ „, ,,. (^-1) 
^2 < ■ ^o\xaXj\/A 

\\S2{Ts + S,)S2\ 



< 1/2 

Then for e > and r oc eL-^/"^, there exists a TCP map E, consisting of L unitaries 
Uj on n + 1 qubits of the form 

U, = cxp [-tT,{Hs + m t^'^ +Ts® + S,)) 

and single qubit measurements, such that for any state \F) e Si, 

Tr[(|0)(0| ® ln+iM\F i){Fi\)] = 1 - 0(e) 
as e — ^ 0"*" , with total simulation time 



i5/2 X L 

T = 



y— 



^A\xo\) j-[ \xj 
Proof: 

The proof is constructive. E can be described by a loop of L cooling steps, labeled 
by the index j (starting at j = L). At the beginning of each iteration, the bath is 
measured in the logical basis. If it is in state |t), then the transition to the ground 
state has already occured, so wc effectively terminate by decrementing j ^ j — 1 
and continuing to the next iteration. Otherwise, wc evolve under Hamiltonian Hj = 

Hj + V + Sj, where Hj = Hs + w^^ls 8) |t)(t |s, V = Ts ® a^, and 6j represents an 
error term in the simulation. The bath energy Ug'^ (defined explicitly below) is near 
uij. The time evolved under this Hamiltonian is Tj = + 0{r^)), and the 

associated unitary evolution is labeled Uj . We then decrement J — > J — 1 and continue 
to the next iteration. 

Written in pseudocode, E may be summarized as 
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For j= L, L-1, ... 1 

Measure bath qubit 

If bath is down: 
Apply U_j 
Measure bath 

Step j is associated with the following trace preserving, completely positive (TCP) 
map: 

E,{p) ^ uMUpimiu] 
+lt>(tlHt)(tl 

where represents the operator (I5 (81 and likewise for tt)(tl- We define 

M = El o E2... o El. Given the above assumptions, the lemmas below prove that the 
algorithm works as expected. The proofs of the lemmas are included at the end of 
this section. 

Lemma 7.2 (Fidelity of the cooling step). Let Uj = exTp{—iTjHj) be the unitary 
evolution associated with cooling state Assume that (7.1) holds. There exists a 

time Tj = 20o|aoa | + Oir^)) and w^' e {ivj - A/4,(.jj + A/4) such that 

U,\ji) = \0t)+Oir) 

where the bound 0{r) is uniform over all j . 

Lemma 7.3 (Preservation of the lower bands). Define 

A'/ji_i=Span{|A:;)U=o 
Then under the assumptions of Lemma 7.2, for any state \v) e 

Uj\v)^Mf^iU,\v) + ^-0{r) 

Lemma 7.4 (Projective mapping of TCP map). Define 

Mj = A//©Span({|Ot)}) 

Then for MjpMj = p, 

Ejif) = XjPj-i +Rj 

where Pj-i is a density matrix satisfying Mj_ipj_iAfj„i = Pj-i, Aj = 1 — Tr[i?j], 
Tr[i?,] = . 0{r), and \\R,\\ ^ ^ ■ 0(r). 

Lemma 7.5 (Success of the Deterministic Algorithm). Given the assumptions 
in the previous lemmas, let p = |F^)(i^^| for some \F) e Si, and assume that 
r = Oo/A (X Then 

Tr[A/oE(/7)] = 1 - 0(e) 
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as e — > 0. 



Since Mq = |0)(0| ® In+i, Lemma 7.5 proves the first claim in Theorem 7.1. 
As shown in the proof of Lemma 7.2, unitary Uj requires an evolution time Tj ^ 
where 51 = fio|2:oa;j|(l + 0{r)). The simulated time required for the algorithm is 
therefore 

L 

L 



2Vt^\x(^ \ ^ \xj 

ye^\x,\)2-.^\x,\ 
where we used the fact that fig ^ = (rA)"i = A^^i^/^ • 0{t^^). 



□ 



7.1. Reduced Algorithm:. It is possible that the decomposition \ G) = Xj \ j) 
contains values of Xj that are exponentially small in n, meaning that the simulation 
time Tj = 0(exp(n)). We therefore discuss conditions under which it is valid to ne- 
glect the cooling of such states, leading to an improved run time. In doing so we 
derive bounds on the allowed error in the preparation of the fiducial state 4). 

Suppose that we choose to skip the cooling of states |j) such that \xj\ < i], for 
some parameter r/. Defining S'l = Span{|j) : \xj\ > 77} ® H^, we may write 

\Fi) = S[\Fi)+f^\F^i) 

Given p — \F ^){F i|, wc may compute 

p^il-lh^F' i){F' l\+R^ (7.3) 

where \F' I) G S[, and R± is an error term with rank 2 and norm ||-R_l|| < 

Defining as the complement of S'l in S'l , we can write ^2 = ^2 + . From 
(7.1) and the triangle inequality, 

\\S[iV + 6)S',\\ ^ 0{Ar) 

\\S',{S)S',\\^A{l/2 + 0{r)) 

As these are the only assumptions necessary to prove Lemma 7.3, for sufficiently 
small r it still holds for the reduced space (M^Lj^)' = ^^/-i riS'i- Combining this with 
Lemma 7.2, we see (as in its proof) that Lemma 7.4 also holds with respect to the space 
Mj = Mj n S[ . We may then apply Lemma 7.5 to the density matrix p' = |F' 1 1 
and get the same fidelity e for the reduced algorithm E' = Ei-^ o ... o Ei^,, where 
ii < 12... < 2l' enumerate the eigenstates in S'l. 

Since E' is applied to p and not p' , we sec that the reduced algorithm fidelity is 

Tr[MoE'(p)] = (1 - |/^nTr[A/oE'(p')] + Tr[A'/oE'(i?i)] 
= (l-|/^n(l-0(e)) + Tr[MoE'(i?^)] 
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Because E' is a composition of L' single qubit measurements and unitary evolutions, 
as seen in the proof of Lemma 7.5 we have that ||E'(i?_L)|| < L' ■ ||i?_L|| = L' ■ 0(|/j^|), 
where L' = d\m{S'i)/2 — 1 < L. Furthermore, as Mq is a rank 2 projector. MqE,'{R±) 
has rank at most 2, so we conclude that 

|Tr[MoE'(iZi)]| < rank(MoE'(fli)) • p/oE'(i?i)|| = L' ■ 0(|/i|) 

Using this result in (7.5), as long as — 0{e/L'), we may still achieve fidelity 
1 — 0{e). Notice this allows us to treat any error in the preparation of \F l) as 
contributing to 

Since we are only cooling states such that \xj\ > rj, the timing of the algorithm 
is then 



< O 

= O 



L' 



Ve77A|a;o| 



In the case where |F) = \G) we see that = X]|a;j|<r) I^jI^ - ^ ' l/j-l - 

L^/'^rj and in order to maintain an infidelity 0(e) it is sufficient to set rj = t/L^/"^. 
This gives 



r = o 



e^Alxol 



7.2. Proofs of Lemmas 7.2-7.5. The proof of Lemma 7.2 is the most involved. 
First, we analyze Hj in the subspace 5'i Mb, and show that it has eigenstates near 
\j I) ± |0 t), with energies near cub ± floxoxj. In order to do this we compute the self 
energy operator 'E,p(z) for the manifold of energies near tUj, and show how the value 

(7) 

of ujg may be calculated to account for energy shifts associated with V. We then 
account for the static error term J, and the possibility that V may couple between 
Si and 82- 

Proof of Lemma 7. 2: 

We begin by analyzing the case where 5 — and V = SiVSi = ilo\G){G\ (g) (Jx- 
Note that Si = (Span{|fc)}^^Q) (g) Hs is then invariant under Hj = Hs + tXt 
|(n+i) _j_ which will simplify our analysis. In order to understand the dynamics 

of Hj, we compute the level shift operator (6.1) to find an effective Hamiltonian for the 
eigenstates of Hj with energy near coj. Let P be the eigenspace of energies between 
A_ = ojj — A/2 and A+ = ojj + A/2. Given that the eigenstates of Hs in ^2 are at 
least A away from those in Si, and joj^' — Wj| < A/4, we have 

P = \ji){ji\ + \Ot){Ot\cSi 

g = ^ |fc t)(fc 1 1 + E 1^ ^X^^ ^ I + ^2Q52 

k^O k^3 

PH,p = ^^^\QW^\ 

l/-r!o|G)(G|®a, 
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Using (6.1), we may now compute 

Ep(z) = PHjP + PVP + PV — ^ — TT^yP 

' ' z- Q{Hj + V)Q (7.6) 

= PHjP + PVP + PV (Gq(z) + Gq{z)VGq{z) + ...) VP 

where Gq{z) ~ z-qh Q ' Since the projector into Si commutes with P,Q,Hj, and 
V, we may replace all operators by their projections into 5'i above. This simplifies 
Sp(z), giving 



Ep(z) = c.,m)0-;| + 4^'Vt)(otl 

+ QoP \G) (a, + noTx (G| Gq{z) \F) a,) a^z) {F\ P 
where cr* is the bath operator 

a,(z) = (1 + (r!o (G| Gq{z) \G) a.,f + (^io (G| Gq{z) \G) a^f + ...) 
As a bath operator, one has that 

(G|GQ(z)|G)=.g,(z)U)(;|+.go(;^)|t)(tl 



9j 



{z) = {Gi\GQ{z)\Gi)^Y. 



\Xk 



|2 



5o(^) = (Gt|GQ(z)|Gt>=E' " 



so 



therefore 



(17o (G| Gq |G) a^f' = {^Ig^goT (I iXi I + I tXt I) 



(7.7) 



"■'^'^ = 1 o2 ^ w I + I t)(t I) (7.i 

Since P\G i) = Xj \j |), P |G t) = xq |0 t), we may use (7.7) and (7.8) to get 

Ep(z) = c.,b-;)(j;| + ^(^')|ot>(otl 



2 p |G;)(Gtl + |Gt)(G;| 

1 - "o5j.9o ^ 

- ^ogj\G t)(G t I + f^o5o|G ;)(G ; l) P ^^'^^ 

\ ^^0 / \xj\'^VLQgo{z) xqXj \ 

w^lV l-f^o5j(^)3o(^) V 2;oa:, |.to|21]o5j(2) / 

where the above matrices are written in the {\i 1) , |0 t)} basis. 

Using the above expression, we define the effective Hamiltonian -ffe// = ^p('^b^)- 
We identify the diagonal elements of the second matrix above as the level shift, and 
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observe these are composed of even powers of fig- If the diagonal terms in i?e// were 
equal, it would cause coherent oscillations \j I) o |0t) at a rate 2fl w 2rio|a;oa;j |. 
We may find such that this is the case, by solving the degree L + 1 polynomial 
equation 

((1 - nlg,{z)go{z)){z ~ lj,) + nl{\xo\^g,{z) - \xj\^go{z))) 

Effectively, we are adjusting the value of so that the even order energy shifts 
induced by V are canceled out. Note that we need a solution for w^' that is contained 
in {uij — A/A,u}j + A/4). Using the fact that fio < A/4, as well as \gj\, jgol < 3^ 

for z ~ cli^' G [ojj — A/4:,u}j + A/4], it is not difficult to show that the left hand 
side is negative for uj^^ ~ ujj — A/4 and positive for w^'' = tUj + A/4. Since the 
left hand side is smooth over this range, by the Intermediate Value Theorem a root 
exists within {ujj — A/4, + A/4). If our computed root is <5a;^^ off from the exact 
solution, from (7.9) we see the two states retain a splitting O(Jcjg^). The necessary 
accuracy (Jw^' can thus be incorporated into the error term SiSSi, as long as Sujg^ = 
0{\\SiSSi\\) = no\xox j\-0{r). 



We therefore assume that w^-' has been chosen so that the diagonal terms in (7.9) 



are equal at z = ujg^ , which means 

H,ff = ^*s^ + n{\ji){ot\ + \ot){ji\) 

where 

n= ^"^f^ =noxox,{l + 0{r^)) 



Lu*^ = + »oixox,- r°/^^^r°^^" = + n ■ o{t) 

and we assume that |a;o/a;j| < 0{1). 

Heff has eigenstates \v±) ~ (|j \.) ± |0 t)) with energy ujg ± $7, so dynamics 
under H^ff for time tj = ^ would map state \j |) — > |0t)- To show that evolution 
under Hj achieves the same mapping, we use Theorem 6.1 and Corollary 6.6 to show 
it has eigenstates and energies near those of H^ff- To do this, we must determine an 
error bound for ||I]p(z) — Heff\\. 

As in (7.6) above, since H^ff ~ T,p{LUg^) we have 

Ep(z) - H^ff = PV (Gq{z) - Gq(l.1^^)) VP (7.10) 

where 

for |0s) e Q. Since jw^'' — Wjj < A/4, the eigenvalues of QHjQ are at least 3A/4 
away from ijjj. As \\V\ \ < A/4, by Lemma 6.2 we then conclude that Es G {— 00,10 j~ 
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A/2]U[a;j+A/2,oo) for all s. Gq{z) is therefore analytic for z £ [u!j~A/A,L0j + A/4\, 
so we may compute its Taylor series expansion about oj^"* G (ujj — A/4,LUj + A/4): 



G, 



For some Zg G [uij — A/4, + A/4], between uj^g and z. Since |zs — E's] > A/4, we 
conclude that 



||GQ(z)-GQ«)ll<k 



<r- U _ 



(7.11) 



||Ep(.)-F.,,||<|z-4^)|(ffl)' 
= |z-c.l^'|(4r)2 

As above, the spectrum of Hef / is contained in [c, d] , where c = w^' — 51 (1 + 0(?')) , 
d ~ uj^'g + (1 + 0(r)). In Theorem 6.1, we consider only values of z in [c — 7, d + 7] 
(7 is the error in the eigenvalues of H^f f , compared to Hj ) . Thus we can determine 
7 self-consistcntly by solving 

j=\z-Ji'^\{Arr 

for z = d + 7 and 2 = — 7. To leading order in r, this gives 

j^n- 0{r) (7.12) 

(in fact 7 = SI • 0{r^), but the following result holds for (7.12) as well). Applying 
Theorem 6.1, we have that the two eigenvalues of Hpff, E± = ± fi, are 7 close to 
the eigenvalues of Hj + V. The relative error in the energy difference (£'+ — E^) is 
therefore 0(7/^!) = 0(r). 

Now Corollary 6.6 can be used to show that the eigenvectors of -ffe// arc close to 
the corresponding eigenvectors of Hj . In the notation of that corollary, we can define 
P' = \v+){v^\, and see that v = 0,1] = 2il. Denoting the analogous eigenvectors and 
eigenvalues of Hj by \v±) and E±, using (7.12) we see that 

\{v+\v+)f^{i+\P'\v+) 

>('-mi(- 

= l-0(r2) 

and likewise for {v^\v-). From this and Theorem 1 we conclude that 

\v±) = \v±) +0{t)^ i=(b- ;) ± |0 t)) + 0{r) 

E+-E- = {E+ - E^) (1 + 0{r)) = 20(1 + 0{r)) 

For time evolution Tj = ^^-e ~ C''--'^^) iniplics the statement of Lemma 7.2. 
To complete the proof, we account for the case when V 7^ SiVSi or (5 7^ by including 
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the effect of these terms in ||I]p(z) — H^ff\\. As long as (7.f 2) still holds, we conclude 
that (7.13) is still valid. The full Hamiltonian is now Hj = Hj + SiVSi +5eff, where 
Seff accounts for the terms we previously neglected. Specifically, 

Seff^5 + iV-SiVSi) 

We wish to compute the bound ||S]p(z) — Heff\\, where now Sp(z) is defined with 
respect to the perturbation SiVSi + Seff (see (7.14) below). As before, we have 
Heff = Sp(cj^^)|_5_^^^^p, with I]p(z)|^^^^^o defined as in (7.6). Suppose \\^p{z) - 
Ep(z)|^ j.^.=oll ~ ^' ■ triangle inequality, 

||Sp(z) - HeffW < y + ||S]p(z)|^^^^^o - HeffW 

We could then repeat the previous analysis to compute 7, and get 7 = Q,-0{r) + 0{'^'). 
The results of Theorem 6.1 and Corollary 6.6 could then still be applied to get (7.13), 
as long as 7' also satisfies (7.12). Below we show that this is the case, as long as (7.1) 
is true. 

Including 5eff in (6.1), we see that 

I]p(z) = PHjP + PVP + PSeffP 

+ P{SiVSi + 5eff) . {SiVSi + <5e//)?'^^^ 

z-Q{H,+S^VS^+5eff)Q 

In order to bound all terms proportional to Seff, we use the relation {A — B)~^ = 
A-^ + A-^B{A - B)-^ to get 

^ GQ{z) + GQ{z)Seff- ^ 



z-Q{Hj + SiVSi + Seff)Q 'z~Q{Hj + SiVSi + Seff)Q 

where Gq{z) = z-q{h +SiV Si)Q • This allows us to write: 

Sp(z) - 5]p(z)b^^^^„ = PSeffP + PSeff ^ SeffP 

" z-Q(Hj + SiVSi+Seff)Q 

+ PSeff ^ —SiVSiP + h.c. (7.15) 

'z-Q{H,+SiVSi+Seff)Q ^ ' 

+ PSiVSiGgSjeff ,^^SiVS,P 

" z~Q{Hj + SiVSi + Seff)Q 

We now bound this difference. The operator Q{z — Hj — SiVSi)Q can be di- 
agonalized in blocks of Si and ^2. As before, for z e [ojj — A/4:,u!j + A/4], within 
both Si and ^2 this operator has eigenvalues with magnitude at least A/2. In the 
notation of Lemma 6.5, we may define A = Q{z — Hj — SiV Si)Q, B ~ QSeffQ, 
G, = A/2,G2 = A/2, so that {A S)- = ^_q,,^^s?vS.^s.,,)q - ^^^^^^ = 

\\S,QSeffQS,\\ll^ < ||SA//S,||/A, i?x - \\SiQSeffQS2\\/A < ||Si^e//S'2||/A, by 
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Lemma 6.5 one may show that 

Salh 0(i?x) 

z-Q{H,+SiVSi+5eff)Q 

l-^i — . ^^..^ , , .^ ^ilh 0(1) 

-S2|hO(l) 

Writing aU terms of (7.15) in Si, S2 blocks, we have for z e [ujj — A/4, + A/4] 

PSeff^A-{ 0{Ri) 0{Ry,) ) 



Q 


z - Q{H, ^ 


h^iF^i - 

Q 


F Seff)Q 


z - QiHj - 


hSiVSi - 


F Seff)Q 




Q 




z - Q{Hj - 


hSiVSi - 


^Seff)Q 



z - Q{Hj + SiVSi + <5e//)g " V 0(i?x) 0(1) 

SiVSiP = A • ^[j' ' 

P5iy5iGQSi<5e// = A • ( 0(ri?i) 0(ri?,) ) 

With these components, using (7.15) one may calculate 7' = ||Sp(z)— I]p(z)|^ //=oll ~ 
0(A(i?i +RI,)). We need ||Sp(z) - Sp(z)|^_^^^^g|| = n ■ 0{r), so we require 

Ri = 0(r • n/A) 
Rl = 0(r • n/A) 
R2 < 1/2 

where the last inequality comes from the bound on S^fj necessary to use Lemma 6.5. 
One may check that these statements are satisfied by (7.1). 

□ 

Proof of Lemma 7.3: 

As before, we have Hj = Hs + Hb, Hj = H + V + Sj. Define Pq as the eigenspace 
of H with energies ujq^uji, ...ojj. This corresponds to the space ^f-i ^ Si mentioned 
in the lemma. In the language of Corollary 6.7, it corresponds to A_fc = Afc+ = LOk 
and A as defined for Hs- The proof comes in two steps. We define the intermediate 
Hamiltonian H'^ = H + SiVSi, with an eigenspace Pi corresponding to energies 
within {ujQ - A/8, ujq + A/8) U {uji - A/8,wi + A/8) U ... U (w^ - A/8, ujj + A/8). 
Likewise, P2 is the eigenspace of Hj of energies within {ojq — A/4, ujq + A/4) U (wi — 
A/4, wi + A/4) U ... U {ujj - A/4, Wj- + A/4). The proof follows by showing that (up 
to an error 0(r^)), any state in Pq is in Pi, and any state in Pi is in P2. This will 
imply that a state in Pq undergoing evolution Uj will remain in Pq. For simplicity of 
notation, for all equations below let \vi) represent a normalized state in P;. 

Since H'^ and Hj are block diagonal in ■S'l and ^2, as long as. Pi C Si it is 

sufficient to reduce our analysis to Si. This holds if S2{V + S)S2 does not change 
the energy of S2 states by more than A/2, as implied by Lemma 6.2 and (7.1). 

Considering only Si, (7.1) implies the bound H^i {h'^ - Si\\/ A = 0{r). By 

Corollary 6.7, we have that 



{vQ\Pi\vQ)^l~j-0{r^) 
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Writing |iio) = a \vi) + b \v^'^ where Pi \v^'^ = 0, one can easily show that 



Pibo) = v/i-.rO(r2)|«i) 

Equations (7.1) also imply that Pi is energetically separate from Qi = 1 — Pi by 
at least A' = A/4, so that \\Hj - H'j\\/A' = 0{r), and as above, 



P2\vi) = v/1-.r 0(r2)|z;2) 
We can combine these statements to get 

{V0\P2\V0) > {vo\PlP2Pl\vo) 

^il-j-0{r^)){vi\P2\vi) 

Finally, writing \vq) = a \v2) + b |i'2"), the above statement implies 

P2W)^{l-J-0{r^))\v2) 
and by an identical analysis, for any \v2) G P27 there exists \vq) G Pq such that 

P^\v2) = {l-]-0{r^))\v^) 

Notice that the diagonals of Pq are at least as large as those of P2PqP2- Since P2 
is an eigenspace of i/j , it is clear that P2Uj ~ UjP2- Using these facts and the above 
equalities, we compute the bound: 

(fol f/jPoC/, \vo) > {vol U]P2PoP2Uj \vo) 
= {vo\P2U]PoUjP2\vo) 
= il-J-0{r^)f{v2\Po\v2) 
= il-j-Oir')f^l-j-Oir')-f-Oir^) 
= l-]-0{r^) 

where in the last line we use the fact that r cx L^^/"^. Since Qo = 1 ^ Po we get 

\WAvo)\? <3-0{r^) 

Writing Uj \vq) ~ PqUj |uo) +QoUj |wo), we conclude the proof noting that Pq = ^if^i 
and that the bound 0{r^) is dependent only on the ratio VIq/ l^. 

□ 

Proof of Lemma l.J^: 

Since the operation Ej starts with a bath measurement and since the Mj projector 
commutes with the bath projectors ||)(| | and |t)(t I, we may assume without loss of 
generality that 

p = A//M^/+Po|Ot)(Otl 
where = Span{|fc ^Yk=Q- Furthermore, since 

t>(0 t I) = |0 t)(0 t I = M,-i£^,(|0 t)(0 t 
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by the linearity of TCP maps it suffices to analyze the component of p within Afj*. 
We may therefore assume that p = M'-'pM'^. Since p is a density matrix, we have 
that 

P = ^Pi\vi I I 
I 

where \vi \) S and I is a sum over at most j + 1 = dim(Mj^) terms. Each \vi 1} 
may be decomposed into components parallel and orthogonal to \ j \.): 

\vi i) = ai \ j I) +bi \ v^) 

where G ^^f-i- By Lemma 7.3, we have that 

Uj\vt)=Mf_,U,\vt) + ^i-0{r) 

Likewise by Lemma 7.2, 

t^.|ji> = |0t)+O(r) 

= M,_i[/, |j ;) + 0{r) 

Since M^_^ C Mj^i, we conclude that 

U,\vi i) = M,_i(7, i) + v^-0(r) 
Finally, since Ej is a linear operator, we see that 

= Y,Pi {U,\vi i){vU\U]) 
I 

= M,_,E,{p)Mj_,+Rj 

where Rj is the sum of all terms proportional to 0{r). Using the triangle inequality 
and the fact that the pi sum to 1, we see that \\Rj\\ = V7 ' 0{r). Since Rj is a 
sum of at most {j + 1) operators, each of rank 2, we see that rank(i?j) < 2{j + 1). 
Therefore \Tr[Rj]\ < rank(Aj) • \\Rj\\ = f^^ ■ 0{r). Since a projection of a density 
matrix is proportional to a density matrix, we may write Mj-iEj(p)Mj^i ~ AjPj_i, 
with = TiiAIj^iEj{p)Mj_i] = Tr[Ej{p) ~ Rj] = 1 - //z . o{r). 

□ 

Proof of Lemma 7. 5: 

E is defined by the chain of TCP maps, 

¥. = Eio E2... o El 
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The initial state of the system and bath is described by the density matrix Pl = 4- 
){F II, where \F I) E AIi^. Repeated apphcation of Lemma 7.4 gives 



= El o E2... o El-1 [XlPl-1 + Rl 

— Eio i?2... o i?L_2 {^lXl-iPl-2 + ^lRl-i + El^i{Rl) 

= (^Y[ p° + 

where MqPqMo = pQ, Rtot represents all other terms, and Tr[i?tot] = 1 — (^nfe=i 
since E is trace-preserving. 

We will bound the infidelity, 1 - Ty[MqE{pl)] = 1 - (llLi ^fe) - T^AMoRtot], by 
showing that 1 — Tr[MoRtot] are small. First, consider the quantity 

y = log (rife^i ^kj = J2k=i log(A/c)- By Lemma 7.4, we see that Afc = 1 - Tr{Rk)- 
Given that | log(l - x)| < 2\x\ for \x\ < 1/2, we conclude that for |Tr(7?fc)| < 1/2, 

L 

|2/|<^|log(l-Tr[i?fc]) 

k=l 
L 

<5^2|Tr[i?,]| 

k=l 
L 

k=l 



. 0(r) 



Thus, to leading order in r, 
1 - 



n ^fe ) = 1 - e« = . 0(r) 



To show the second term is small, we must bound Rtot, which is the sum of all 
error terms: 

k=l \j = k+l J 

From the simple form of Ej (see (7.2)), we see that 

k 

El o ... o Ek-i{Rk) - ^J^kA] 

where 

A, = itxtK^.imi) • {u,+i\m\) ■ ■ ■ {Uk-i\m\) 
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for 2 < j < k, and 



Ai = {Ui\m\) ■ mmi) ■ - ■ iUk-i\m\) 



Since Aj is a product of projectors and unitaries, we must have that ||74ji?fcAj|| < 
so by the triangle inequality and Lemma 7.4 it follows that 

\\E, o ... o Ek-iiRk)\\ < k\\Rk\\ = . 0{r) =^ 



\\Rt,t\\=J2k'/'-0{r)=L'/^-0ir) 
fc=i 

Finally, we note that since Mq is a projector of rank two, Moi?tot^^o also has 
rank at most 2. By the cyclic property of the trace, we conclude that 

\TT[MoRtot]\ - \TT[MoRtotMo]\ < mnk{MoRtotMo) ■ \\MoRtotMo\\ - L^^^ ■ 0{r) 

Combining the two results, we have that 

l-Tr[AfoE(pi)] =i5/2.o(^) 

Thus, as long as r = 0(e/L^/^), the algorithm succeeds with infidelity 0(e). 

□ 

8. Extension Analysis. We now discuss an augmentation of the previous cool- 
ing technique which docs not require knowledge of the overlaps Xk describing the 
fiducial state, \F l). It is described in detail in the article, though we summarize it 
here. We assume that the system Hamiltonian Hs has the form 

Hs = PiHsPi + P2HSP2 

where Pi is the eigcnspace of Hs with energy between (uji — Su) and {ui + Slu), |0) 
is the nondegenerate groundstates Hs with energy wq = 0, and P2 is a projector into 
the space orthogonal to Span{|0)}®Pi. To relate to notation in the previous section, 
we define the projectors Si = (|0)(0| + Pi) ^ 1, S2 = P2 ^ 1- Since these act trivially 
on the bath, as a slight abuse of notation we will sometimes refer to Si and S2 as 
operating on the system Hilbert space alone. We define the spectral gap between 
|0),Pi andPs: 

A = rRm^uJi,E,\E-uji\ : E e Spcc{Hs\sj] (8.1) 
In full, the unperturbed Hamiltonian is 

H^Hs®{\C){C\ + \R){R\-\L){L\) 
+ u:ils<»{\R){R\ + \L){L\) 

where |C) , \R) and \L) are orthogonal basis vectors for the bath Hilbert space. We 
start by preparing a fiducial state. 



\FC) = h\lC)+U\F^C) 
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where |1) G Pi, (1| F±) = 0, and we are given a lower bound for The algorithm 
proceeds by simulating the evolution of Hamiltonians H + X , where X satisfies 

X = Ts®{\C){B\ + \B){C\) 

\B) = ^i\L) + \R)) (8.3) 

n\l)=PiTs\0) 

where by phase convenstion n is real. Again, although 1 1 1 1 = fio is a known quantity, 
we are only given a lower bound H* for f2. 

The algorithm is probabilistic, and involves a single evolution step for time t ~ ^, 
followed by a measurement of the bath. If the bath is measured in state \B), then 
the desired transition |1C) — > \0 B) could have occured. We verify this by applying 
the bath unitary \B) o \L) ,\D) ■(-> |_R), then evolving under H + Y for time r = 
where 

Y = nols®i\L){R\ + \R){L\) (8.4) 

A measurement of a bath transition \L) — > \R) would indicate that the system is in its 
ground state, while in all other cases a transition \ipL) — > \ipR) is suppressed by energy 
conservation. If either the first or second bath measurements fail, we reinitialize the 
system and start again. 

As before, we again show that given some bounds, the algorithm is robust against 
simulation errors and coupling between Si and 5'2: 

r = rio/A < 1/8 
Suj <r- fl{V) 
|(0|T5|0)i2 



UJl 



< r ■ Vl{V) 



\Si5Si\\ myi (8.5) 

A A 

\Siiy + 5)S2\\'' n{v) 

< r ■ 



A2 ^ ■ A 

\S2{V + S)S2\ 



A 



< A/2 



where n{V) = fta for V ^Y, and n{V) ^ {l\Ts |0) = 1^ for F = X. As seen below, 
for the probabilistic scheme to succeed with fidelity 1 — 0{e), we must scale r as 
0{\fi\e'/'). 

Lemma 8.1 (Fidelity of the Unitary Evolutions). Let U{t) = e-*^(«+^+'5) and 
assume (8.5). Then 

U{t) \\C)= cos(0t) |1 C) -ism{(j)t) |0 B) + 0{r) (8.6) 

where (j)t = tQ[1 + 0{r)). The error term in (j)t and in (8.6) is uniform over r. 
Likewise, let U„{t) = £-'^(^+'^+'5). Then 

C/„(t)|OL) = cos((/)i,) |0 L) -isin(0„)|Oi?) + 0(r) (8.7) 

where 6,, — rrJn. 
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Leivima 8.2 (Verification Step). 



max(7/;L| [/] (I5 ® \R){R\) \ip L) = ©(r^) 



(8 



,8) 



where the maximum is taken over all normalized system states \ip) such that {ip\0) — 0. 

Theorem 8.3 (Success of the Probabilistic Scheme). Say that r = 0(|/i|e^/^) 
as e — > 0^, and that the Ut simulation time t is sampled randomly within the range 
[j^,^], where fl* < fl. Then the verification step accepts with probability p„ = 



state is Psuccess = 1 ^ 0{e). Since < rA, the average simulation time (T) satisfies 



The proof of Lemma 8.1 is analogous to the proof of Lemma 7.2. We first analyze 
the success of the unitary evolutions under the assumption the most unwanted terms 
are zero, and show that it leads to the desired outcome. We then bound the effect of 
the unwanted terms on the unitary evolution. 
Proof of Lemma 8. 1 : 

We begin by proving the first statement of the lemma, for V = X = Ts ® 
(\B){C\ + \C){B\). As in the proof of Lemma 7.2, instead of analyzing H + V + 6 we 
start by looking at the evolution oi H + SiV Si, then obtain a bound on the errors 
caused by jg// — V ~ SiVSi + 6. Define P as the eigenspace of H with energy 



in [wi - A/4,wi + A/4]. Notice that the projector P is P = Pi «) \C){C\ + |0)(0| (g) 
i\D){D\ + \B){B\), where \B) = ^(|L) + |P)), \D) = ^J\L) - \R)), and that P C Si. 



Before we calculate the self energy operator S(z) at z = wi, we note the following 
relations: 



P = |0)(0| (|L)(L| + |P)(P|) + Pi ® \C){C\ 
= \0}{0\^{\D){D\ + \B){B\)+Pi(^\C){C\ 
Q = Pi® {\L){L\ + |P)(P|) + |0)(0| ® \C){C\ + S2 
PHP = a;i|0)(0| ® {\D){D\ + |P)(P|) + Pi^s^'i ® \C){C\ 
QHQ = (PiHsPi + S2HSS2) ® i-\L){L\ + \R){R\) 

+wi(Pi + S2) » {\L){L\ + \R){R\) + S2HSS2 <E} \C){C\ 
PVP - P [T5 ® {\C){B\ + \B){C\)] P 
^n{\lC){OB\ + \OB){lC\) 



The next term required in (6.1) is the unperturbed Green's function, Gq{z) — ^qJquq 
Since P d Si and S'lV^S'i, H are both block diagonal in Si and ^2, we may ignore 
the S'2 component of Gq{z): 





SiGq{z)Si 



Pi ® (|L)(L| + |P)(P|) 



[z - Loi) - PiHsPi ® (|P)(P| - \L){L\) 
+ i|0)(0|®|C)(C| 
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so that 

5iGQ(a.i)5i = ® {\L){L\ \R){R\) 

+-|0)(0|®|C)(C| 

Notice that {B\ SiGq{oji)Si \B) = 0. From the definition of SiVSi = S'lTsS'i® 
{\C){B\ + \B){C\), we immediately observe that 

SiVSiGq{uji)SiVSi = —SiTs\OmTsSi ® 

Multiplication by Gq{uji) again produces a term proportional to \D). Since V\D) =0, 
this implies that the series (6.1) with perturbation SiVSi truncates at second order 
in V. Therefore i?e// = may be computed to all orders as 

Heff = PHP + PVP + PVSiGq{loi)SiVP 

= t^i|0)(0| ® {\D){D\ + \B){B\) + PiHsPi ® \C){G\ 
+n{\lG){QB\ + \OB){lC\) 
, l(0|Ts|0)p 



-mm 



The system Hamiltonian is written Hs — PiHsPi + S2HsS2, where the spectrum 
of Hs\pi is contained in {uJi — Suj , uji + 6uj) . The state |1) oc PiTs |0) is not necessarily 
an eigenstate of Hs, but = (1| Hs |1) is contained in (wi — 6uj, uji + Suj). We write 
the projector into the remainder of Pi as P{ = Pi — |1)(1|. To see that H^ff produces 
the desired evolution, we rewrite it as 

Heff - {\OB){OB\ + \1G){IC\) + iuJi\OD){OD\ + P[HsP[ ® \G){C\) 

+ n{\lC){OB\ + \OB){lC\) ^gg^ 

+ (K - ^i)ii>(ii + piHsimi + h.c.) ® \G){c\ + MIiBl\oB)m\ 

LOl 

Observe that if we neglect the terms on the third line of (8.9), Heff has eigen- 
vectors \v±) = ^ (lie) ± \QB)) with eigenvalues ± f], which exactly produce the 

desired evolution (8.6). Since ||P{(i?s ~ uJiP[)P[\\ < \\Pi{Hs -a;iPi)Pi|| < 5uj, by 
Lemma 6.2 all other eigenvalues of the approximate H^ff are in (wi — 5uj,uji + 5uj). 
The eigenvalues wi ± are therefore non-degenerate and energetically separated from 
the rest of the spectrum by a gap VI — Suj. This fact will allow us to use Theorem 6.3 
below to show that, up to an error of order 0{r), \v±) correspond to eigenvectors of 

Heff- 

The terms in the third line of (8.9) are bounded by ft ■ 0{r). To see this, 
note that ^''"^"^''^^^^^ < 51 • 0{r) is already an explicit assumption. The bound for 
((wi* - wi)|l)(lT+ PiHs\l){l\ + h.c.) = Pi{Hs - LOiPi)Pi - Pi{H, - LOiP[)P[ comes 
from the fact that (i?s - ujiP'i)P[\\ < \\PiiHs - ujiPi)Pi\\ < 6uj ^ fl ■ 0{r). By 
invoking Lemma 6.2 and Theorem 6.3 we conclude that H^ff has eigenvectors |±) 
with eigenvalues loi ± fl ■ (1 + 0{r)) such that | (±| ±)p > 1 — O(r^), and that the 
rest of the spectrum of H^ff is ■ {1 + 0{r)) away from these eigenvalues. 

The rest of the proof of (8.6) is now identical to the argument in Lemma 7.2. Using 
the bound for fip = l|Ts||i in th^ case when Seff = we bound ||S(z)|^ — H^ffW 
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for 2 G [wi — 51 • (1 + 0{r)),LL!i + ■ {1 + 0{r))] using the Taylor's expansion of 
— — -- — — E?— — , . Then, usina; Lemma 6.5 we bound the error in Efz) obtained 

by neglecting 6eff ~V — S\VS\ + 5, and show that it is equal to • 0{r) under our 
assumed bounds (8.5). Since ||S(z) — i?e//|| is still sufficiently small, we conclude by 
Theorem 6.1 and corollary 6.6 that + F + j has eigenvalues wi ± $7 • (1 + ©(r)), 
and that these eigenvalues correspond to (|1 C) + |0i3)) + 0[r). 

The proof of the second statement is nearly identical to the first. Noting that 
now V = 17ols ® + wc have PVP ^ no\OR){OL\ + h.c. and PVQ = 0. 

Assuming that S = 0, the level shift operator S(^) is now exactly equal to PHP + 
PVP, so Heff ~ S(a;i) satisfies 

Heff = Loi (|OL)(OL| + \QR){QR\) + PiHsPi ® \C){C\ 
+ 17o(|0L)(0i?| + |Oi?)(OL|) 

which clearly has eigenvalues uji ± 51o corresponding to :^(|0i) ± |0i?)), with the rest 
of its spectrum in {uji — 5lj,lji + 5uj). H^ff therefore produces the desired evolution. 
The rest of the proof, in which wc bound — i?e//||, again continues in the same 

way as in Lemma 7.2, with the substitution of fig in place of fl. 

□ 

Proof of Lemma 8. 2: 

The proof of this analogous to the proof of Lemma 7.3. As before, let P represent the 
eigenspace of H with energy contained in (— oo,a;i — A]. Notice that P corresponds 
only to bath states in state \C) or \L). Q = 1 — P corresponds to the eigenspace of 
energies within [wi,c»). In the language of Theorem 6.3, we have A_ = — c», A+ = 
wi — A/2, and A defined as in (8.1). 

Given (8.5), by the triangle inequality we conclude that ||y+(5||/A = 0{r). Define 
P as the eigenspace of H=H + Y + 6 with energy in (A_ , A+) = (— oo, cji — A/2). By 
Theorem 6.3, for any state \v) e P, we see that 




This implies that P\v) = -^/l — O(r^) \v), where \v) represents an arbitrary normal- 
ized vector in P. Likewise, for \v) S P, {v\ P\v) = 1 — 0{r^). Since P is an eigenspace 
of ff, for Uy = e^p{-iTH), PUy = UyP. Noting that {^\ PPP [V') < (V'l P IV') for all 
we conclude that for any \v) G P, 

{v\ UlPUy \v) > {v\ UlPPPUy \v) 

= {(v\PU}\p(UyP\v) 



= {l-0{r^)){i\P\i) 

= {I - 0{r^)f = \ - Oir^) 

Let |'0 L) be given such that ('010) = 0. By examinig the spectrum of H , one sees 
that L) G P, and that the eigenspace Ms ® Span{|i?)} is contained within Q, so 
that the operator I5 ® |^)(^| < Q = (1 ~ P)- From the above inequality we conclude 
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that 

L\ Ul (Is ® \R){R\) IV L) < (V L\ Ul (1 - P) \i' L) 

= l-{iPL\UlPU,\i:L) 

By an identical argument, we may prove the second statement of the claim as well. 

□ 

Proof of Theorem 8.3: 

Given initial state \FC) and time evolution C/t(r) (as defined in Lemma 8.f), the 
probability of a verification event is given by 

p„ = {FC\ ulx^XUt \FC) 

where 

X = \R){R\U,\L){B\ 

and Uy is evaluated for time Likewise, the probability of the system being in the 
ground state after the verification has occurred is 

Pys ^ {FC\ujx^\0){0\XUt\FC) 

where |0) is the ground state of Hs- Along with finding p^, we wish to calculate the 
success probability of the algorithm conditional on a verification event, p^^y = p^s/Pv 
These three probabilities are functions of the parameter (j>, where (f>/T = + 0{r)) 
is half the energy splitting of the eigenstates \v±) ~ ^(|1C) ± \0B)) used in the 
evolution Utir). 

By the first result of Lemma 8.1, for any time evolution Utir) we may write 

Ut{T) \FC) = h (cos(0t) \IC) - zsin(0O |0B) + 0{r)) + f^Ut \f^C) 

Since Lemma 8.1 implies coherent oscillations between |1C) and |0i?), it must be 
that (0B| Utir) |/_lC) = 0(r), so that (0B| C/t(r) \FC) = -i{fi sm{(j)t) + 0{r)). We 
mention that the bound 0{r) is independent of r, i.e. as r — > 0+ there is a constant 
c> such that I (OB I Utir) \FC) + i{fi sin(0t)l < c • r for all r. 
Likewise, the second result of Lemma 8.1 implies that 

^"(^=2^) I0^) = |0^)+O(0 

so 

X\OB) ^ \0R) + O{r) 

For system states such that (-01 0) = 0, by Lemma 8.2 we have that 

max (0L| Ul\R){R\Uv \iJj L) = 0{r^) 
{■4>\i>)=i 

and since U^\R){R\Uv is a projector, it must be that Ul\R){R\Uy {t/j L) — 0{r), where 
the error bound is uniform over all \tp). We conclude that for all composite states 
\iPsb) such that (-0ss| OB) ^ 0, 

X \^sb) = Uy {Ul\R){R\Uy) ■ i\L){B\^PsB)) = 0(r) 



36 



Algorithmic Cooling of a Quantum Simulator 



and that this bound is the same over all such IV'ss)- Since {OB\Ut{T)\FC) = 
— sin((/)t) + 0(?')), we may write 

Utir) \FC) = -lifi sin(0O + 0(r)) |0 B) + c(t) IVsbW) 

where |c(t)| < 1 and (V'ss(t)I OB) = for all r. 
Applying X gives 

X[/t(r) \FC) = sin(</)0 + 0(r)) |0i?) + 0(r) 

-t{h sin(0t) + |0i?) + B{cbt) \^±{t) R) 

where 0) = for all r. As the bounds derived from Lemma 8.1 and 9 are 

independent of r, we have that A = max = 0{r) and B = max |_B((/)t)| = 0{r). 

We may now directly compute p„ and pys ■ 

pMt) = l/i sin(0O + + mcl^t)? 

Pvs{(bt) = |/isin(0t)+^(0t)l' 

so the success probability for a given 0t is 

. . ^ |/isin(0O+A(^t)P 
^^1''^'^*^ |/isin(0O+^(0t)P + |S('/'t)P 

Suppose that when sampling over values of 0t, with probability at least 1 — e/2 
we have Ps\v{4>t) > (1 ~ e/2). Then Psuccess > (1 — > 1 — e for < e < 1, which 

is the desired result. In terms of the relation above, this condition equivalent to 



\hsm{<j>t) + Aicl,t)\ > |i?(0t)|y^-l 
As l/i sin(0t) + A{(j)t)\ > l/i sm{(j)t)\ — A and |-B(0t)| < B, this relation is satisfied if 



A + B 



|s.„(*)l>^^^O^^J (8.10) 

Hence we obtain an infidelity at most e as long (8.10) is violated with probability 
at most e/2. Note that if (pt is sampled uniformly over a range larger than tt/2 (as 
ensured by our sampling scheme for r), the probability that |sin(0()| < c for some 
small number c > is pfaii = 0(c) as c — )- 0. We require pfaii < e/2 in (8.10), which 
is satisfied for r = Od/ije^/^). This gives success the bound Psuccess > 1 — e stated in 
the Lemma. Using the same argument, we see that to have |/i sin((/)t)| > 2A ^ 0{r) 
with probability at least 1/2, we only require r = O (|/i|), so under the more stringent 
scaling we may also conclude that the verification probability py{(j)t) is greater than 

l/ip/4 with probability 0(1). This gives the desired scaling, paccept = O (j/Tp) ■ 

□ 

9. Universality. The following is derived from results in [36, 38, 35]. Using the 
notation of [38], |0") represents a state of n qubits, each initialized in the qubit state 
|0). The universality of both QSC schemes follows immediately from this claim, and 
the fact that 1 and 2-local unitaries are universal [33] . 
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Theorem 9.1 (Universality of QSC). Let U = UlUl-i---Ui be composed of L 
one and two-qubit gates on n qubits. There exists a 5-local Hamiltonian Hg on n + L 
qubits, whose ground state \ri) tracks the history of the unitary evolution: 

L 

=^C/zC/z_i...[/i|0")|l'o^-') 

Furthermore, Hg has a subspace Si that contains \r]) , composed of L+1 nondegenerate 
eigenstates. These states are resolved by at least A, with A^^ = O (riL^^. The state 
\G) = |0"0^) G Si has spectral decomposition 



L 



where \k) is the kth excited state in Si. Finally, Si is also invariant under the single 
qubit operator Tg ~ i^oln ® |0)(0| ® li_(„+i), which satisfies SiTgSi = f2o|G)(G|. 

Using the results of Theorem 7.1, Theorem 9.1 implies that we may produce the 
history state of U using total simulation time T = ^ ■ O ( j^^ v^i°s(^) ^ ^ g^j^^^ ^^^^^ 

of IliJsIlT = O — ^'°g(^) ^ , If we coneatenate ^-j^L identity operations to the 

definition of U, we see that \ri) is then e-close (with respect to the trace- norm) to the 
state U |0"). We conclude that any computational problem on n qubits that may be 
solved at poly(n) cost using standard quantum computation may also be solved by 
QSC at a cost of poly(n). 

The ground state of Hs can also be produced with the alternative scheme. As 
seen in the proof of Theorem 9.1, in the language of that scheme we may define 
|0) = I?/), |1) = \k = 1) as any other eigenstate in the low energy subspace, and 
\F) ~ \G) = |0"0^). In this case, the probabilistic scheme produces \r]) at an average 

cost of \\Hs\\{T) = O (-^^)- We could produce the state U |0"0^) by adding M 

identities to U, then measuring if the clock states one of |l-^0^+^) through |l^^)- 
Since the scheme is already probabilistic, the cost of producing (710") would only 
change by a constant multiple factor. 

Proof of Theorem 9.1: 

For this we use Kitaev's clock Hamiltonian [36], which acts on a system of n + L 
qubits. The first n qubits represent the actual computation, and start in the fiducial 
state |0"). The L other qubits are the 'clock' qubits (denoted by '^) that keep track of 
the evolution of the system. In the notation of [35] , assuming a characteristic energy 
scale uj, the Hamiltonian is 

where 

^f™p«t = ^|l)(l|,:®|0)(0|5 (9.3) 

1=1 

ensures that the initial clock state |0^)^ is associated with the initial computational 
state |0"), 

L-l 

H^iock - E |01X01|m+i (9-4) 

(=1 
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gives an energy cost for not being a valid clock state |l'0^~')'^, and 



Hprop 2 



1 = 1 

where for 1 < / < i 

H, = l^ (|100)(100|?_i,,+i + |110)(110|?_i,,+i) 

- C/i|110)(100|^_i,;,+i - C//|100)(110|^_i,,;+i 

correspond to the tracked evolution of the computational bits. Hi and are sim- 
ilarly defined, but with clock qubits and L + I omitted. In the above notation, 
the subscripts refer to action on a specific qubit, and imply that other qubits are left 
unchanged by the operator. 

Notice that Hinput, Hprop and Hdock are each positive semidefinite. Define Siegai 
as the null space of Hdock, which is compose of states the form |-) \lY = \-) |l'0^^') '. 
Since Hdock commutes with Hprop and Hinput, we see that eigenstates of Hs in 
the complement of Siegai have energy at least 2w. We now reduce our analysis to 
Siegai, since our desired subspace 5'i will be contained in Siegai and will describe 
energies less than 2cj — A. As in [38], within Siegai we apply the change of basis 
W = J2iLo UiUi-i...Ui Hprop and Hinput are then mapped to 



n 

= Eiixii.®io)(or 

i=l 

^p™p = ^i^El'^-ix^-ir- 



z^i 

-\i~iw-m-i\ 



Since it is tridiagonal, it is not hard to show that Hprop has eigenstates of the form 




N> = N> I \l V-^ > COS ( il±lIlpL\ 



1 

with eigenvalue E'^, = 1 — cos (^-^fij for < k < L. We see that these eigenvalues are 
separated by at least Eiq = Ei — Eq = 2sin^(j^), which is greater than ^^^_^iyi for 
L > I. Furthermore, the largest eigenvalue El is bounded by 1 — cos (^x^pr) ; ^o Siegai 
is separated in energy from its complement by at least lj ■ -q^^, thereby justifying 
our reduction to Siegai- 

Within Siegai, for each logical state |a;) on the n computational qubits, define the 
space by its projector 5* = |a;)(a;| (g) J2^Lo lOd'^- '^x commutes with Hinput and 
Hprop, so each forms an invariant subspace of HsW. Furthermore, within we 
see that 

SlH.nputSl = N{x)\x){x\ ® |0)(0r 
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where from the previous equation 



and N(x) is the number of I's in the binary expression for x. 
Within S'*, the Hamiltonian is then 

SlW^HsWSl = SlW^LuHprop + AiH,nput)WSl 

= \x){x\^ (^^Y^EklvM'' + Nix)A,\0m'^^ 

Scahng Ai/w as h ■ 2n(L+i)'^ ^ ' ^w/n, we sec that ||Aiifi„put|| < h ■ ujEiq, where 
ujEio is the minimum eigenvalue spacing of Hprop- Hence within each space S* we may 
treat Ai_ff,j„p„t as a perturbation to ujHpj-op, which for small h is well approximated 
by first order perturbation theory: 

Ek^ Ek + Ai {vk\ H^nput \vk) (1 + 0{h)) 

= Ek + At,—-^ cos^ — (1 + 0(h)) 

* ^ L + l \2{L+l) J ^ ^ " 

Since N{0) = 0, we see then that the eigenstates of H input +E[ prop in Sq^ are separated 
in energy from the other invariant subspaces by at least Ai-^-ijCOS^ ( 2(^+1) ) ~^ 
0[h)). We conclude that the eigenstates in are spectrally resolved by the gap A, 
with 

A-i = A7i(L + l)cos-2 (^^^—\ {l + 0{h)) 



\2{L+l) 



- ■ O inL"^) 



With this information, we may define S\ by its projector 

Si = wsi^w'' 

Si is spanned by the eigenstates 
\k)^W\On\r,k) 



L 



'fY^cos(^^]u,U,^i...Uim\l) 



L 



1=0 



with energy Wfc = w • (1 — cos y^fij), where |?/) is the fc = state. Furthermore, these 

eigenstates are non-degenerate, and gapped by A as defined above. By the definition 
of the eigenstates and (9.6), we conclude that the state \G) = |0") |0^)'' = |0") \l = 0)'' 
is contained in Si, and that it satisfies (9.1). We note that 



so the operator T5 = |0)(0|5; leaves the space Si invariant, with SiTsSi = |G')(G'|. 
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